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Abstract 

Most supersymmetric models predict new particles within the reach of the next 
generation of colliders. For an understanding of the model structure and the 
mechanism(s) of electroweak symmetry breaking, it is important to know the 
masses of the new particles precisely. The measurement of the mass of the scalar 
partner of the top quark (stop) at an e + e~ collider is studied. A relatively light 
stop is motivated by attempts to explain electroweak baryogenesis and can play 
an important role in dark matter annihilation. A method is presented which 
makes use of cross-section measurements near the pair-production threshold as 
well as at higher center-of-mass energies. It is shown that this method does not 
only increase the statistical precision, but also reduces the influence of system- 
atic uncertainties, which can be important. Numerical results are presented, 
based on a realistic event simulation, for two signal selection strategies: us- 
ing conventional selection cuts, and using an Iterative Discriminant Analysis 
(IDA). While the analysis of stops is particularly challenging due to the possi- 
bility of stop hadronization and fragmentation, the general procedure could be 
applied to many precision mass measurements. 
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Abstract. Most supersymmetric models predict new particles within the reach of the next genera- 
tion of colliders. For an understanding of the model structure and the mechanism(s) of electroweak 
symmetry breaking, it is important to know the masses of the new particles precisely. The mea- 
surement of the mass of the scalar partner of the top quark (stop) at an e + e~ collider is studied. A 
relatively light stop is motivated by attempts to explain electroweak baryogenesis and can play an 
important role in dark matter annihilation. A method is presented which makes use of cross-section 
measurements near the pair-production threshold as well as at higher center-of-mass energies. It is 
shown that this method does not only increase the statistical precision, but also reduces the influ- 
ence of systematic uncertainties, which can be important. Numerical results are presented, based 
on a realistic event simulation, for two signal selection strategies: using conventional selection cuts, 
and using an Iterative Discriminant Analysis (IDA). While the analysis of stops is particularly 
challenging due to the possibility of stop hadronization and fragmentation, the general procedure 
could be applied to many precision mass measurements. 

PACS. 14.80.Ly Supersymmetric partners of known particles - 95.35.+d Dark matter (stellar, 
interstellar, galactic, and cosmological) 



1 Introduction 

Supersymmetric particles are likely to be produced 
and observed in high-energy proton-proton collisions 
at the LHC. However, it will be difficult to confirm 
their identity as superpartners of the known Standard 
Model particles and to measure their properties pre- 
cisely. For this, one needs experiments at a linear e + e _ 
collider such as the proposed ILC at y/s = 500 GeV. 

An experiment at the ILC will be able to make 
many precise measurements from which particle prop- 
erties, and ultimately, the outlines of a particle physics 
model may be inferred. Due to the high statistical pre- 
cision expected at the ILC, the optimization of the 
systematic errors is of particular importance. We have 
studied one specific example, the extraction of the 
mass of a scalar top quark from cross-section mea- 
surements near threshold. We have devised a method 
which reduces most systematic uncertainties and leads 
to a potentially very accurate measurement of the stop 
quark mass. This method is general and could be ap- 
plied to other particles that are pair-produced in an 
e + e~ collider. 

The method relies on the comparison of produc- 
tion rates at two different center-of-mass energies, and 
knowledge of how the cross-section varies as a function 
of i/s and the mass of the particle. 
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We have chosen the case of a light scalar top with 
a mass not much higher than the mass of the light- 
est neutralino since production of this particle was al- 
ready extensively studied in an ILC context. It was 
concluded that a conventional approach to the mea- 
surement of the stop quark mass culminated in an un- 
certainty of about 1 GeV pQ. The new method im- 
proves substantially on this result. The presented re- 
sults are preliminary and being finalized [SJ. 

For this analysis, we have performed realistic sim- 
ulations of the signal and backgrounds, and used two 
techniques to separate the signal from the background. 
The first technique is based on conventional sequential 
cuts, while the second employs an Iterative Discrimi- 
nant Analysis (IDA). Furthermore, the hadronization 
followed by fragmentation of the stop has been in- 
cluded and we have carefully studied the systematic 
uncertainties arising from this and other sources. 

There are theoretical motivations for studying a 
light stop quark with a small mass difference. Specif- 
ically, we evoke a scenario within the Minimal Super- 
symmetric extension of the Standard Model (MSSM) 
which is able to explain the dark matter density of the 
universe as well as the baryon asymmetry through the 
mechanism of electroweak baryogenesis [I] . 

A small mass difference between the stop and the 
lightest neutralino can help to bring the dark mat- 
ter relic density into the observed region due to co- 
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annihilation between the stop and the neutralino. For 
this mechanism to be effective, the typical mass dif- 
ference is rather small, Am — mz — m^a < 30 GcV. 

' L l Al 

The dominant decay mode of the stop is ti — ► c\\, 
resulting in a final state with two soft charm jets and 
missing energy. 

Previous methods to determine the scalar top quark 
mass were discussed for the SPS-5 benchmark (mi = 
220.7 GeV) [3] and results are summarized in Table [1] 
For the cosmology motivated benchmark with m~ ti = 
122.5 GeV and m^o = 107.2 GeV, an experimental 
precision of Ami = ±1.0 GeV was obtained pQ, and 
about ±1.2 GcV including theoretical uncertainties. 
The following study investigates the same signal sce- 
nario and it is based on the same background reactions 
and event preselection. 

Table 1. Comparison of precision for scalar top mass de- 
termination for the SPS-5 benchmark = 220.7 GeV). 



Method 


Am- U (GeV) 


luminosity 


Polarization 


0.57 


2 x 500 ft -1 


Threshold scan 


1.2 


300 m- 1 


End point 


1.7 


500 ft -1 


Minimum mass 


1.5 


500 m- 1 



2 Mass Determination Method 

This method proposes to derive the stop mass from 
measurements at two center-of-mass energies, one mea- 
suring the stop production cross-section near the thresh- 
old (th), and the other measuring it at a center-of-mass 
energy where the cross-section has approximately a 
peak (pk). Using both measurements leads to a can- 
cellation of systematic uncertainties in the mass deter- 
mination. A parameter Y is defined as 

Y _ ^th — -6 t h _ Cth £th ^ th 
V p k — -Bpk fpk £pk £pk 

where TV is the total number of expected events af- 
ter event selection and B the number of correspond- 
ing background events, a is the stop production cross- 
section, e the selection efficiency, and C the luminosity. 
The center-of-mass energies 260 and 500 GeV have 
been chosen. Near the threshold, the production cross- 
section is very sensitive to the stop mass. 

In this study we assume that the ILC will operate 
primarily at \fs — 500 GeV with a total luminosity of 
C = 500 fb _1 , and a small luminosity of C = 50 fb -1 
will be collected at y/s = 260 GcV. Table[2]summarizes 
the expected production cross-sections. The detector 
response was modeled with the SIMDET package in- 
cluding the LCFI vertex detector concept. 

The relation of the observable Y and the stop mass 
is used to determine the stop mass with precision. For 
example a variation of Y by 3% in a realistic scenario 
would lead to an uncertainty of the stop mass Ami = 
0.2 GeV as illustrated in Fig. [TJ 

3 Sequential-cut Analysis 

In order to cancel the systematic uncertainties to a 
large extent with the described method, the same se- 



Table 2. Cross-sections for the stop signal and Standard 
Model background processes for i/s = 260 GeV and yfs = 
500 GeV and different polarization combinations. The sig- 
nal is given for a right-chiral stop of rrif = 122.5 GeV. Neg- 
ative polarization values refer to left-handed polarization 
and positive values to right-handed polarization. 



Process 


a (pb) 


at y/~S = 


260 GcV 


a (pb) 


at y/s = 


500 GcV 


P(e-)/P(e + 


0/0 


-.8/ + . 6 


+ .8/-. 6 


0/0 


-.&/+. 6 


+ .8/-. 6 


hi* 


0.032 


0.017 


0.077 


0.118 


0.072 


0.276 




16.9 


48.6 


1.77 


8.6 


24.5 


0.77 


zz 


1.12 


2.28 


0.99 


0.49 


1.02 


0.44 




1.73 


3.04 


0.50 


6.14 


10.6 


1.82 


eeZ 


5.1 


6.0 


4.3 


7.5 


8.5 


6.2 


qq, q ^ t 


49.5 


92.7 


53.1 


13.1 


25.4 


14.9 


ti 


0.0 


0.0 


0.0 


0.55 


1.13 


0.50 


2-photon 
p T > 5 GcV 


786 






936 







0.11 



0.1 




0.09 H — .— i — .— H — 1 '' 1 i 1 — I— ■ — 1 

121.5 122 122.5 123 123.5 

stop mass (GeV) 

Fig. 1. Example of mass uncertainty derivation from the 
uncertainty of the observable Y. 

quential cuts are applied for the y/s = 260 and 500 GeV 
analyses. Details of the event selection are given in Ta- 
ble [3] and the results are given in Table |4j 



Table 3. Selection cuts for y/s = 260 and 500 GeV. Also 
listed are the selection efficiencies optimized for right-chiral 
stop quarks. 



Variable 


= 260 GeV 


V7 = 500 GeV 


number of charged tracks 
visible energy E vis 


5 < "tracks < 2 5 
0.1 < -Eyis/v/a < 0-3 


5 < "tracks < 20 
0.1 < -Eyis/v/s < 0.3 


event long, momentum 


Pi/Ptotl < 0.85 


[Pt/Ptotl < °- 8 5 


event transv, momentum pf 


15 < p T < 45 GcV 


22 < p T < 50 GeV 


thrust T 


0.77 < T < 0.97 


0.55 < T < 0.90 


Number of jets iVj eta 


"jets > 2 


"jets > 2 


extra-jet veto 


E jot < 25 GcV 


E jot < 25 GcV 


charm tagging likelihood P c 


Pc > 0.6 


P c > 0.6 


di-jet invariant mass mjj 


tti 2 . < 5500 GeV 2 or 

JJ 2 2 
m.. > 8000 GcV 2 

11 


Tn 2 . < 5500 GeV 2 or 

J i 2 

m 2 > 10000 GcV 2 


signal efficiency 


0.340 


0.212 



Table 4. Numbers of generated events, and expected 
events for the sequential-cut analysis at i/s = 260 and 
500 GeV for total luminosities of 50 ft)- 1 and 200 ft- 1 
with unpolarized and polarized beams. 





y/5 = 


260 GcV 






500 GeV 






generated 


L = 50 


fb- 1 




£ = 200 


fb" 1 


P(e-)/P(e+) 




0/0 


8/-. 6 




0/0 


.8/-. 6 





50,000 


544 


1309 


50,000 


5170 


12093 


w+ w 


180,000 


38 


4 


210, 000 


16 


2 


zz 


30,000 


8 


7 


30, 000 


36 


32 


Weu 


210,000 


208 


60 


210, 000 


7416 


2198 


eeZ 


210,000 


2 


2 


210, 000 


< 7 


< 6 


qq, q ^ t 


350,000 


42 


45 


350, 000 


15 


17 


tt 










180, 000 


7 


7 


2-photon 


1.6 X 10 6 


53 


53 


8.5 X 10 8 


12 


12 


total background 




351 


171 




7509 


2274 
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4 Iterative Discriminant Analysis (IDA) 5 Systematic Uncertainties 



The Iterative Discriminant Analysis (IDA) is applied 
to increase the discriminant power between signal and 
background compared to the sequential-cut-based anal- 
ysis, and thus reduce the statistical uncertainty in the 
stop mass measurement. Figures [2] and [3] give the re- 
sults of expected number of background events as a 
function of the signal efficiency. The chosen working 
points have efficiencies of 38.7% and 41.6% for the 
y/s = 260 and 500 GeV analyses, respectively. Table [5] 
lists the corresponding expected background events. 
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Fig. 2. IDA: Expected background events as a function of 
the signal efficiency for C — 50 fb _1 at y/s = 260 GeV. 




0.3 0.4 0.5 
Efficiency 

Fig. 3. IDA: Expected background events as a function of 
the signal efficiency for C = 200 fb _1 at y/s = 500 GeV. 



Table 5. Numbers of generated events, and expected 
events for the IDA at y/s = 260 and 500 GeV for total 
luminosities of 50 fb _1 and 200 fb~ x with unpolarized and 
polarized beams. 



■/s = 260 GeV 



■/J = 500 GoV 



£ = 200 fb" 



P (e~)/P(e^) 



.8/-. 6 



0/0 



50.000 



50.000 



1V T W" 

ZZ 

Wev 

eeZ 

qq> q =t ■ 
tt 

2-photor 



180,000 
30,000 
210,000 
210,000 
350,000 



210, 000 
30, 000 
210, 000 
210, 000 
350, 000 
180, 000 
5.5 X 10 6 



20 
1719 
< 7 
18 
1 



510 

< 6 



Both the sequential-cut-based analysis and the IDA 
method lead to a small statistical uncertainty resulting 
in Ami < 0.2 GeV and thus systematic uncertainties 
are particularly important to evaluate. Three classes 
of systematic uncertainties are distinguished: 

— instrumental uncertainties related to the detector 
and accelerator: detector calibration (energy scale), 
track reconstruction efficiency, charm-quark tag- 
ging efficiency, and integrated luminosity. 

— Monte Carlo modeling uncertainty of the signal: 
charm and stop fragmentation effects. The Peter- 
son fragmentation function was used with e c = 
-0.031±0.011 (OPAL). For e b = -0.0041 ±0.0004 
(OPAL) and e b = -0.0031 ± 0.0006 (ALEPH) an 
average uncertainty of 15% was taken, and a factor 
2 improvement at the ILC has been assumed, lead- 
ing to Aei 1 = 0.6 x 10~ 6 where ej x — eb(mh / m mi ^) 2 . 
Fragmentation effects and gluon radiation increase 
the number of jets significantly and the importance 
of c-quark tagging is stressed in order to resolve the 
combinatorics. 

— neutralino mass 108.2 ± 0.3 GeV. 

— theoretical uncertainties on the signal and back- 
ground. Some improvement compared to the cur- 
rent loop calculation techniques is assumed, and an 
even larger reduction of this uncertainty is antici- 
pated before the start of the ILC operation. 

Tables [6] and [7] list the systematic uncertainties for 
the sequential-cut analysis and the IDA. The system- 
atic uncertainty using the IDA method from detec- 
tor calibration (energy scale) is larger. This is because 
the sequential-cut analysis pays particular attention 
to cancellation of this uncertainty between the two 
analyses at the different center-of-mass energies. 

Table 6. Sequential-cut analysis experimental systematic 
uncertainties on the signal efficiency. The first column in- 
dicates the variable that is cut on. The second column con- 
tains the expected systematic uncertainty for this variable 
based on experience from LEP. The third column shows by 
how much the signal efficiency for y/s = 260 GeV varies as 
a result of varying the cut value by this uncertainty. The 
fourth column shows the same for y/s = 500 GeV, and the 
fifth column lists the resulting error on the observable Y. 







rcl. shift 


on signal eff. 






260 GeV 


500 GeV error on Y 


energy scale 


1% 


3.7% 


3.1% 0.6% 


^tracks 


0.5% 




negligible 


charm tagging 


0.5% 




taken to be 0.5% 


luminosity 




0.4% 


0.2% 0.4% 


charm fragmentation 


0.011 


0.3% 


0.8% 0.6% 


stop fragmentation 


0.6 X 10~ 6 


0.6% 


0.2% 0.7% 




0.3 GeV 


3.8% 


3.0% 0.8% 


background estimate 




0.8% 


0.1% 0.8% 


Table 7. IDA 


experimental systematic uncertainties. 






rcl. shift 


on signal eff. 




variable 


260 GeV 


500 GoV error on Y 


energy scale 


1% 


3.4% 


1.3% 2.3% 


^tracks 


0.5% 




negligible 


charm tagging 


0.5% 




taken to be 0.5% 


luminosity 




0.4% 


0.2% 0.4% 


charm fragmentation 


0.011 


0.1% 


0.6% 0.5% 


stop fragmentation 


0.6 X 10~ 6 


0.1% 


0.8% 0.7% 




0.3 GeV 


3.7% 


1.6% 2.2% 


background estimate 




0.3% 


0.2% 0.1% 
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6 Mass Determination 

The assessment of the achievable stop mass precision 
is based on the statistical and systematic uncertainties 
on the observable Y (eq. (1)) as summarized in Table[5] 
The IDA method has a smaller statistical uncertainty, 
and also a smaller background uncertainty due to a 
smaller number of expected background events. The 
expected stop mass uncertainty is inferred from the 
uncertainty on Y as given in Table [9] 

Table 8. Summary of statistical and systematic uncertain- 
ties on the observable Y. 



for Y 



sequential cuts 



statistical 
detector effects 
charm fragmcntatk 
stop fragmentation 

background estima 



3.1% 
0.9% 
0.6% 
0.7% 
0.8% 
0.8% 



2.7% 
2.4% 
0.5% 
0.7% 
2.2% 
0.1% 



of experimental systc 
of experimental error 



1.7% 
3.5% 



3.4% 
4.3% 



theory for signal ( 



Table 9. Estimated measurement errors (in GeV) on the 
stop quark mass. 



error category 



sequential cuts 



ti 



statistical 
sum of expo 
beam spectr 
sum of expe, 
sum of all C3 



cntal systema 
and calibrate 

and th. errors 



0.19 
0.10 



0.24 
0.42 



0.17 

0.21 
0.1 
0.28 
0.44 



7 Cold Dark Matter (CDM) Interpretation 

The chosen benchmark parameters are compatible with 
the mechanism of electroweak baryogenesis pQ. They 
correspond to a value for the dark matter relic abun- 
dance within the WMAP bounds, /2cdm^ 2 = 0.109. 
The relic dark matter density has been computed as in 
Ref. pQ . The assumed benchmark parameters changed 
slighty (larger slepton masses assumed) and ^cdm/i 2 
changed from 0.1122 PQ to 0.109. In the investigated 
scenario, the stop and lightest neutralino masses are 
= 122.5 GeV and m^o = 107.2 GeV, and the stop 
mixing angle is cos6*j = 0.0105, i.e. the light stop is 
almost completely right-chiral. The improvement com- 
pared to Ref. pQ regarding the CDM precision determi- 
nation is shown in Fig.Uand summarized in Table. fTTJl 



0.14- 




123 



Fig. 4. Computation of dark matter relic abundance 
OcdmIi taking into account estimated experimental errors 
for stop, chargino, neutralino sector measurements at the 
future ILC. The black dots correspond to a scan over the la 
(A\ 2 < 1) region including the total expected experimen- 
tal uncertainties (detector and simulation) , the grey-dotted 
region includes also the theory uncertainty, and the light 
grey-dotted area are the previous results PQ. The red star 
indicates the best-fit point. The horizontal shaded bands 
show the la and 2a constraints on the relic density mea- 
sured by WMAP. 

8 Conclusions 

Scalar top quarks could be studied with precision at 
a future International Linear Collider (ILC) . The sim- 
ulations for small stop-neutralino mass difference are 
motivated by cosmology. The precision mass determi- 
nation at the future ILC is possible with a method 
using two center-of-mass energies, e.g. y/s = 260 and 
500 GeV. This method can also be applied to other 
analyses to improve the mass resolution in searches 
for new particles. The precision of two independent 
analysis methods, one with a sequential-cuts and the 
other with an Interative Discriminant Analysis (IDA) 
lead to very similar results. The new proposed method 
increases the mass precision by about a factor of three 
due to the error cancellation using two center-of-mass 
energies with one near the production threshold. In- 
cluding experimental and theoretical uncertainties, the 
mass of a 122.5 GeV scalar top could be determined 
with a precision of 0.42 GeV. The interpretation of this 
benchmark scenario leads to a uncertainty on ficDM^ 2 
of —0.013 and +0.015 which is about a factor two 
better compared to previous results, and comparable 
to current cosmological (WMAP) measurement uncer- 
tainties. With the new stop mass determination, the 
stop mass uncertainty is no longer the dominant un- 
certainty in the X?cdm^- 2 calculation. 



Table 10. Estimated precision for the determination of 
stop mass and dark matter relic density for different as- 
sumptions about the systematic errors. 





Am~ tl (GeV) 




exp. and th. errors 


0.42 


o inq +0 015 

U.1U9_ 013 


stat. and exp. errors only 






sequential-cut analysis 


0.24 


109+ 0012 


IDA 


0.28 


109+ 0012 
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